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Abstract. We present a 39 h integration with the Odin satellite on the ground-state 118.75 GHz line of O2 
towards the region of strongest molecular emission in the Small Magellanic Cloud. Our 3<r upper limit to the 
, O2 integrated intensity of < 0.049 K km s _1 in a 9' (160 pc) diameter beam corresponds to an upper limit on 

the O2/H2 abundance ratio of < 1.3 x 10 -6 . Although a factor of 20 above the best limit on the O2 abundance 
obtained for a Galactic source, our result has interesting implications for understanding oxygen chemistry at sub- 
solar metal abundances. We compare our abundance limit to a variety of astrochemical models and find that, at 
low metallicities, the low O2 abundance is most likely produced by the effects of photo-dissociation on molecular 
cloud structure. Freeze-out of molecules onto dust grains may also be consistent with the observed abundance 
f*^ . limit, although such models have not yet been run at sub-solar initial metallicities. 

^: 

^ \ Key words, galaxies: individual: SMC - ISM: molecules - astrochemistry 

■ 1. Introduction and discuss the implications of our results for astrochem- 
ical models of gas at low metallicities. 

. £^ I Oxygen chemistry in the interstellar medium is still not 

^ '• well understood. Despite having an abundance that is Molecular oxygen has not yet been detected convinc- 
H ; roughly twice that of carbon, oxygen-bearing species are ingly in any source in the Milky Way. Goldsmith et al. 
. - - difficult to observe in the gas phase. Molecular oxygen, ( 2003 report 3<r upper limits on the 2 abundance from 
originally predicted to be almost as abundant as CO un- the SWAS satellite of < 2.6 x 10" 7 in star-forming clouds 
der some conditions (Graedel et al. 11555) 1 . remains elusive, and < 3 x 10~ 6 in cold dark clouds. Pagani et al. 1|2UU3[1 
Although 2 is not a significant reservoir of oxygen in the present improved 2 upper limits from the Odin satellite, 
gas phase (Bergin et al. l2UUU)) . a detection of 2 would be wit h 3<j upper limits to the 2 abundance in 11 sources 
very helpful in distinguishing among the many different rangin g from < 5.2 x 10~ 8 to < 5.7 x 10~ 7 . Goldsmith 
astrochemical models developed to explain its low abun- et al. |2QQ2| report a possible detection in the p Oph A 
dance (see review in Goldsmith et al. 2000). In this Letter, outflow with N(0 2 )/N(H 2 ) ~ 10~ 5 . However, the Odin 
we present upper limits for the 2 abundance towards a upper limit ov er a si milar region of the sky is < 9.3 x 10~ 8 
single line of sight in the Small Magellanic Cloud (SMC) (Pagani et al. I2003fl and so the SWAS detection must be 
considered tentative at best. 

Send offprint requests to: C. Wilson , , , , , , , - c , , . 

* , , . . , „ ,. ,. , , , „. Fxtragalactic sources at moderate redsmtts have their 

Based on observations with Odin, a Swedish- led satellite _ . . . . . ... . . . . 

. , . , , a i . i tv t , . , c -r, , 2 lines shifted away from high opacity regions of the 

project funded jointly by the Swedish National Space Board J " ^ J ° 

(SNSB), the Canadian Space Agency (CSA), the National Earth s atmosphere. Liszt $^SMk obtained an upper 

Technology Agency of Finland (Tekes), and Centre National llmlt to the 2 /CO abundance ratio of 0.1 111 four Seyfert 

d'Etudes Spatiales (CNES). The Swedish Space Corporation galaxies. By searching for 2 in absorption in a system at 

was the industrial prime contractor and is also responsible for z ~ 0. 685 to wards the radio source B0218+357, Combes 

the satellite operation. et al. (1997) obtained a 3er upper limit to the 2 /CO 



2 



Wilson et al.: Searching for O2 in the SMC 




D 50 100 150 

Velocity (km/s) 



Fig. 1. O2 spectrum towards the position SW-1 
(00:46:30.08, -73:21:07.5, J2000) in the SMC. The units 
of the vertical axis are K (T£) and the expected velocity 
of the undetected O2 line is indicated by the horizontal 
bar in the upper right corner. 

abundance ratio of 6 x 10 -3 . If the CO abundance in this 
extragalactic source is similar to that in Galactic molec- 
ular clouds, the 3<r upper limit to the O2/H2 abundance 
ratio would be 6 x 10~ 7 . Thus, even the most sensitive 
extragalactic O2 searches have now been surpassed by the 
Galactic results from Odin and SWAS (Goldsmith et al. 
I20UU1 Pagani et al. lM%|) . 

There have been no previous searches for O2 emission 
from objects with low intrinsic metal abundance. Frayer 
& Brown (1997) present a set of astrochemical evolution 
models aimed at predicting the abundances of key molec- 
ular species, including O2, in high-redshift galaxies. Their 
results suggest that the molecular oxygen abundance is 
enhanced when the metallicity is reduced by even a fac- 
tor of a few relative to the solar abundance. These results 
motivated us to search for O2 emission towards the Small 
Magellanic Cloud (SMC), which has an oxygen abundance 
12 + log(0/H) = 7.96 (Vermeij & van der Hulst I2UU2|) or 
0.16 times the abundance of Orion (Peimbert et al. H993|l . 
At a distance of 60 kpc (Harries et al. I2003f) previous CO 
observations (Rubio et al. I1991f) suggest that molecular 
emission should fill a substantial fraction of the large Odin 
beam (9' at 119 GHz), which makes the SMC the ideal 
target to search for O2 emission from low-metallicity gas. 

2. Observations and Analysis 

We observed the position SW-1 in the SMC (Rubio et 
al. lTWTjl with Odin for 158 orbits (39.2 h on-source) be- 
tween 2004 May 13 and 2004 June 2. The beam diam- 
eter of Odin's 1.1 m telescope is 9' at the 118.750 GHz 
frequency of the O2 li — lo hue with a main beam effi- 
ciency of 0.9 (Frisk et al. I2H0H1 . The observations were 
made in position-switch mode with the off position at 
(+30', 0). The pointing is estimated to be accurate to 15". 
The typical single sideband system temperature of the O2 



HEMT receiver (Frisk et al. 2003) during this period was 
T sys ~ 950 K. The O2 receiver was not phase-locked dur- 
ing these observations, so the true frequency was recovered 
using telluric lines observed when the Odin beam passed 
through the Earth's atmosphere at the beginning and end 
of each orbit (see Larsson et al. 120031 for a detailed de- 
scription of the method applied to the NH3 line) . Satellite 
tests indicate a frequency stability of better than 1 MHz 
(2.6 km s _1 ). The backend was an autocorrelator with a 
spectral resolution of 125 kHz (0.316 km s _1 ) and 700 
channels for a total bandwidth of 87 MHz (220 km s" 1 ). 

The spectrum was fit with a 5th order polynomial base- 
line over the velocity range -35 to 155 km s . A window 
of 20 km s _1 centered on Vlsr = 120 km s -1 , the ob- 
served velocity of the CO J=l-0 line (Rubio et al. [1991), 
was excluded from the fit. No O2 line is detected in the 
final spectrum (Fig.^l, which has an rms noise of 7.3 mK 
(TV) per spectral channel. Using the formula from Pagani 
et al. lj2TH^jl with the CO line width of 15.9 km s" 1 ob- 
served in an 8.8' beam (Rubio et al. 11991(1 . the 3a upper 
limit to the O2 integrated intensity is < 0.049 K km s _1 . 

3. An Upper Limit to the O2/CO Abundance 
Ratio 

We need an estimate of the kinetic temperature to con- 
vert the O2 integrated intensity into an O2 column den- 
sity, N(C>2) (see Pagani et al. 2003 and references therein). 
The Odin O2 beam contains a weak IRAS source (Rubio 
et al.EnSI) and also a 170 /im source seen by the Infrared 
Space Observatory (Wilke et al. 12301 ■ The 60 to 170 /jm 
spectral energy distribution (SED) of this source (n 5 #ll) 
is more peaked towards longer wavelengths than is the 
average SED for the SMC, which Wilke et al. j|2lM|) fit 
with three modified blackbodies of temperature 45 K, 
20.5 K, and 10 K. Thus, a kinetic temperature of 10- 
20 K is probably most appropriate for the region ob- 
served by Odin. A kinetic temperature of T = 20 K gives 
iV(0 2 ) < 2.2 x 10 15 cm~ 2 , while a temperature of 10 K 
gives iV(02) < 1-9 x 10 15 cm -2 . We therefore adopt a 
value N(02) < 2.0 x 10 15 cm~ 2 as the 3<r upper limit to 
the O2 column density in our further analysis. 

At the position SMC-SW1, the CO J=l-0 integrated 
intensity in an 8.8' beam is ~ 1 K km s _1 (Rubio et al. 
I1991|l . Adopting a CO-to-H 2 conversion factor of 1.5 x 10 21 
cm -2 (K km s" 1 )" 1 to account for the lower metallicity 
of the SMC (Wilson 1995) gives an H2 column density 
N(R 2 ) = 1.5 x 10 21 cm" 2 . In the nearby region LIRS36, 
Chin et al. (1998) measured a CO abundance of 1 x 10 -5 . 
Applying this CO abundance to the broader SMC-SW1 
region gives a CO column density JV(CO) = 1.5 x 10 16 
cm -2 . We can also estimate N(CO) directly from the CO 
integrated intensity. Adopting the [ 12 CO]/[ 13 CO] abun- 
dance ratio measured in the solar neighborhood (~ 60, 
Langer & Penzias ll 993(1 and a 12 CO/ 13 CO intensity ratio 
of ~ 12 (Rubio et al. fTMBfl gives a 12 CO optical depth 
of ~ 5. Using the formula from White & Sandell i|1995|> 
and adopting an excitation temperature of 20 K gives a 
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CO column density of 5.3 x 10 15 cm -2 , roughly a factor of 
three smaller than that derived via the CO-to-H2 conver- 
sion factor. For the rest of this paper, we adopt the H 2 and 
CO column densities derived via the CO-to-H 2 conversion 
factor, which correspond to 3er upper limits to the 2 /H 2 
abundance of < 1.3 x 10~ 6 and the 2 /CO abundance of 
< 0.13, both averaged over a 150 pc diameter beam. 

4. Implications for Understanding Oxygen 
Chemistry 

For comparison with chemical models, it is useful to know 
the gas-phase abundances of O and C in the SMC. Dufour 
et al. (TT?M|l derive a value of 12 + log(C/H) = 7.16 in HII 
regions in the SMC. Combining this with the O abun- 
dance from Vermeij & van der Hulst (2002) and in the 
absence of any differential depletion of oxygen relative to 
carbon onto dust grains, the SMC should have a ratio C/O 
= 0.16. This value is significantly smaller than the C/O 
abundance ratio of 0.4 that is commonly assumed as the 
initial gas-phase abundance in astrochemical models. 

4.1. Comparison with Pure Gas-Phase Models 

Goldsmith et al. 12000(1 review the key processes and is- 
sues involved in oxygen chemistry at Galactic metallici- 
ties. Steady-state, pure gas-phase models predict 2 abun- 
dances of 5 — 10 x 10~ 6 , much higher than the current 
observational limits (Pagani et al. I2003f) . One way to de- 
crease the gas phase abundance of 2 is to increase the 
gas phase C/O ratio to near unity via some kind of prefer- 
ential depletion of O onto dust grains, presumably in the 
form of H 2 and other ices. Bergin et al. (2000) discuss 
results for pure gas-phase models with C/O abundances 
ranging from 0.4 to 1. To fit the observed low 2 abun- 
dances, either the gas phase C/O ratio must be close to 
1, or the gas must be "chemically young" , i.e., have spent 
at most 10 5 yr in a high A v environment. This apparent 
youthfulness can be produced by processes such as turbu- 
lent mixing of shielded and unshielded regions (i.e. Chieze 
& Pineau des Forets|l989), which can reduce the 2 abun- 
dance by factors of 100-1000 compared to the steady state 
results. 

The SMC has a C/O ratio of 0.16, much lower than the 
initial Galactic ratio of 0.4. For the Bergin et al. (2000) 
gas-phase model to be consistent with our new limit for 
the 2 abundance in the SMC, either oxygen must be 
preferentially depleted even more strongly (by a factor of 
five relative to the undepleted C/O ratio), or the gas in 
the SMC must be both "chemically young" and suffer pref- 
erential oxygen depletion by a factor of ~ 2.5. Both these 
scenarios seem unlikely, especially since there is still not a 
good justification of why oxygen should be preferentially 
depleted from the gas phase relative to carbon (Roberts 

& HerbstEnni- 

Viti et al. (|2001|l discuss steady-state solutions of gas- 
phase chemical models involving bi-stability solutions. In 
general, if the gas ends up in the high ionization phase, 



then the 2 abundances are low and consistent with cur- 
rent upper limits. In the low ionization phase, the 2 
abundance is too high to be consistent with current up- 
per limits. The bistability region occurs for 0.49 < C/O < 
0.66, with the high ionization phase being the sole solu- 
tion for larger values of C/O. In addition, as the C/O ratio 
decreases, the bi-stability region occurs at lower densities. 
Thus, if the clouds in the SMC have relatively low den- 
sities (riH 2 < 1000 cm~ 3 or lower), the gas could be in a 
steady-state high ionization phase that meets our limits 
on the 2 abundance. 

Spaans & van Dishoeck (2001 ) have run chemical mod- 
els of clumpy clouds where the deeper penetration of the 
UV field affects the chemistry and find 2 abundances be- 
low 10 -7 for 7Jh 2 ~ 10 3 cm -3 . If the gas-phase abundance 
of oxygen is depleted by a factor of two, the 2 abundance 
in their models drops by a further factor of 20. However, 
these models have been run only for solar C/O values. 

4.2. Comparison with Models Including Dust Grains 

Astrochemical models that include the effects of dust 
grains appear promising for explaining the observed low 
abundance of 2 . Roberts & Herbst (2002) use recent lab- 
oratory results for the mobility of H on grains to model 
regions with T = 10 - 20 K and n H2 = 10 4 - 10 5 cm" 3 . 
All their models are consistent with the SMC limits for 
2 abundance at a wide range of evolutionary times. 
However, it is important to note that their models have 
difficulty producing the low Odin limit on 2 in L134N 
(Pagani et al. 127)1131) . 

Viti et al. (|2001|l present a grid of models that include 
the effects of freeze-out of molecules onto grains. These 
models are generally consistent with our SMC upper lim- 
its, whether general freeze-out or selective freeze-out (with 
CO and N 2 desorbing from the grains) is used. However, 
these models use an initial C/O ratio of 0.5, which, for 
the SMC, corresponds to preferential oxygen depletion by 
a factor of 3. Viti et al. ; 20()JJ claim that the selective 
freeze-out models give a chemistry similar to what would 
be obtained with a C/O abundance ratio ~ 1, although 
whether that interpretation is true regardless of the initial 
gas-phase C/O ratio is not addressed in the paper. 

4.3. Comparison with Extragalactic Chemical Models 

Frayer & Brown (JT22Z) combine galactic chemical evolu- 
tion models with gas-phase astrochemical models to pre- 
dict the abundance of 2 relative to CO averaged over an 
entire galaxy as a function of time and metallicity. Their 
models cover a wide range of average metallicity, use pure 
gas-phase astrochemistry, and consider both steady-state 
and "early-time" solutions from a variety of models. They 
also consider selective depletion of oxygen in the gas phase 
by a factor of two, but this does not have a large impact on 
the 2 / CO ratio at low metallicities. Interestingly, at the 
O/C ratio of the SMC (6.3), there is no significant differ- 
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ence in the predicted abundance ratio of O2/CO at early 
times compared to steady state. Both sets of models pre- 
dict O2/CO ratios of 2-3, significantly larger than our new 
observational limit of 2 /CO < 0.13. Thus, even early- 
time gas-phase models appear to have difficulty matching 
our observations of the SMC. 

Frayer & Brown (199TJ) also consider the effects of 
photo-dissociation in low metallicity gas. O2 has a higher 
photo-dissociation rate than CO (van Dishoeck ll988J) and 
so will have a smaller spatial extent inside a molecular 
cloud. For Galactic metallicities and UV fields, Frayer 
& Brown l|1997[l estimate that CO will be self-shielding 
where A v > 1.5 mag and O2 only where A v > 5 mag. At 
lower extinctions, the abundances of CO and O2 depend 
on edge effects in the molecular cloud and are strongly 
non- linear with extinction (see Frayer & Brown (1997) 
for details). If the dust-to-gas ratio scales roughly linearly 
with metallicity (Gordon et al. I2003f) . a molecular cloud 
in a low-metallicity galaxy may not be self-shielding in 
O2 even in its center and may just barely achieve strong 
self-shielding for CO. These effects produce a much lower 
abundance of CO and especially O2 compared to H2 in low 
metallicity gas. Specifically, for the metallicity of the SMC, 
photo-dissociation reduces the O2/CO ratio to a value of 
~ 0.05. Thus, gas-phase astrochemical models which in- 
clude the effects of photo-dissociation at low metallicities 
appear to be consistent with our new observational limit. 

5. Conclusions 

A very deep (39 h) integration on the ground state line 
of O2 with the Odin satellite has failed to detect any 2 
emission towards the region SW1 in the SMC. Our 3er up- 
per limit on the 2 integrated intensity of 0.049 K km 
s _1 translates into an upper limit on the O2 column den- 
sity of N(0 2 ) < 2.0 x 10 15 cm~ 2 . We have combined this 
measurement with published CO J =1-0 data for the same 
region to place a limit on the molecular oxygen abundance 
O2/H2 < 1.3 x 10 -6 and on the O2/CO abundance ratio 
of < 0.13. Although the 2 abundance limit is substan- 
tially higher than the best limit set for a Galactic source 
(< 5 x 10~ 8 , Pagani et al. I2003J) . it is still an interesting 
probe of astrochemistry in low metallicity gas. 

In contrast to the Galaxy, "chemically young" mod- 
els appear to be ruled out for the SMC (Frayer & Brown 
1997). Gas-phase models where the oxygen has been pref- 
erentially depleted from the gas phase compared to car- 
bon do agree with our upper limit for the O2 abundance. 
However, the depletion required for the SMC is a factor 
of two larger than that required to account for O2 limits 
in the Galaxy (Bergin et al. 12QQ0J ), and there is as yet no 
generally accepted mechanism to produce such a depletion 
(Roberts & Herbst EES) • 

Astrochemical models which include the effects of 
photo-dissociation produce an 2 /CO ratio of ~ 0.05 
(Frayer & Brown 11997(1 and are consistent with our 3a 
upper limit. Models that include the effects of dust grains 
also may be consistent with our observed upper limit, al- 



though these models (Viti et al. 120011 Roberts & Herbst 
2002) have not yet been run for reduced metallicity gas. 
If low-metallicity models that included dust grains were 
to fail to produce sufficiently low O2 abundances, then 
the leading explanation for the low O2 abundance in low 
metallicity gas would be the effect of photo-dissociation 
on the structure of the molecular clouds (Frayer & Brown 
I1997fl . This could stimulate renewed attention to the ef- 
fect of clumpy cloud structures and photo-dissociation (i.e. 
Spaans & van Dishoeck 2001 ) in understanding the abun- 
dance of 2 in Galactic sources. 
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